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SUMMARY 


This paper deals with the determination of stress concentration factors 
around a circular hole in a composite laminate. The specific case investigated 
is a four layer (-45°/450/450/-45°) graphite epoxy laminate. The factors are 
determined experimentally by means of electrical resistance strain gages, and 
analytically by using a hybrid finite-element analysis. 


introduction 


In this study, the laminar stress concentrations around a circular hole in 
an angle-ply composite laminate ate determined for the axial tension loading 
case. Of particular interest is the largest value of ctq present at the peri- 
meter of the hole. This study proposes to determine these stresses experimen- 
tally and analytically. For the experimental analysis, electrical resistance 
strain gages are used. The analytic procedure uses the finite-element method 
of a two-dimensional hybrid model with an assumed stress field within the ele- 
ment and assumed displacements at the element interfaces, 


The stress concentration factors around a circular hole in an infinite, 
isotropic sheet have been determined analytically through various approaches 
and confirmed experimentally. For an infinite plate, the stress components 
around the hole are (ref. 1) 
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rhere b is the radius of the hole and CTq is the applied load. The ratio of 
tq/cto along the hole is plotted as shown in figure 1. Obviously, the .maximum 
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of CT 9 is three times a^, and occurs at 9 = + 90°, i.e., at the ends of the 

diameter perpendicular to the direction of tension. 

Due to the increasing use of advanced laminated composites in flight 
structures and other potential applications , the stress concentration around a 
cutout in a fiber-reinforced laminate has been the subject of research by 
several investigators in recent years. Daniel and Rowland (ref. 2) used an 
experimental approach - the Moire technique, and determined the strain (stress) 
concentration around a circular hole in a tension loaded anisotropic plate. 
Hyman et al. (ref. 3) carried some exploratory tests on the same problem. 
Franklin (ref. 4) also investigated the hole stress concentrations in fila- 
mentary structures. By using linear elastic plane stress conditions with the 
help of finite- element method, Ryblcki and Hooper (ref. 5) studied and obtained 
results for boron-epoxy lamina. In a reviewing article (ref. 6 ), Grimes and 
Greimann gave an up-to-date overall picture about the stress concentration 
around a circular hole in a fiber- reinforced composite. Several additional 
references are cited in this article. 

In an experimental study of orthotropic composite materials, Kulkami, 
Rosen, and Zweben (ref . 7) have found that the stress concentration factors are 
a function of the hole diameter, up to a diameter of 2.54 cm (1 in.). They 
observed that the actual number of filaments severed by the hole determined 
the strength of the specimen. 

The present problem of a general angle-ply composite laminate with a 
circular hole is further complicated by the interaction of the individual 
layers. 

EXPERIMENTAL WORK 


Equipment 

The orientation of the strain gages around the holes is shown for each of 
the four specimens in figure 2. The gages were mounted adjacent to the hole 
and were 4.8 mm (3/16 in.) wide, 120 ohm standard foil gages. The specimens 
were mounted in clamp grips and attached to a 90,000 N (20,000 lb.) capacity 
load cell through the use of swivel bearings. The gages were wired into the 
digital strain indicator with the indicator providing three arms of the Wheat- 
stone bridge required in the electrical circuit. The load cell was wired into 
an electrical transducer and calibrated to measure the axial tension applied to 
the specimens. 

The strain gages were applied to the specimens using Eastman 910 adhesive, 
following standard preparation of the surfaces. 


Test Specimen Data 

The four testing specimens were provided by Lockheed-Geolrgia Aircraft 
Company. Their assistance is greatly appreciated. 
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Material: graphite epoxy (Narmco 5209/T300) 

65% graphite fiber, 35% epoxy matrix 
Four layer angle-ply (-45°/45°/45*^-45°) 

Grip tabs of fiberglass epoxy molded integrally with specimens . 
For a unidirectional single layer the macroscopic properties are 


Eqo = 137900 ~ 144795 MN/m^ (20 ~ 21(10)^ psi) 
EgQO = 8274 9653 MN/m^ (1.2 ~ 1.4(10)^ psi) 

G = 4.55 m/r/ (0.66(10)^ psi) 


Specimen No. 

1 

2 

3 

4 

Width 

10.16 cm 
(4 in) 

10.16 cm 
(4 in) 

10.16 cm 
(4 in) 

10.16 cm 
(4 in) 

Hole diameter 

2.5522 cm 
(1.0048 in) 

2.5527 cm 
(1.0050 in) 

2.5527 cm 
(1.0050 in) 

2.5530 cm 
(1.0051in) 


The thickness of the four specimens around the hole was also carefully 
measured. Data were taken at eight stations, the end points of a horizontal 
diameter, a vertical diameter, and two more diameters which bisect the hori- 
zontal and vertical directions. The results are shown in table 1. 

It can therefore be concluded that the assumed thickness 0.6350 ram 
(0.0250 in) is quite reasonable. 


Testing Procedure 

First, the testing specimen was mounted in the upper grips of the loading 
device. The loading indicator and the strain indicator were zeroed and cali- 
brated. Then the other end of the specimen was mounted in the lower grips of 
the loading device. After the specimen was loaded up to the 100% load (2224 N 
or 500 lb), the load was then released. This was repeated six times in order 
to eliminate the strain gage error due to strain hardening. An increiaent of 
20 % of the maximum load was used each time, and the corresponding strain read- 
ing was then taken. The same steps were followed for the other three specimena 


Testing Results 

The data obtained from the strain gage testing was the values for S 0 at 
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four different locations around the hole. These values are given in taole 2, 
and displayed graphically in figure 3. 

The strain gage results can be easily repeated and showed very good 
stability with repeated loadings. The values obtained at 444.8 N (100 lb) of 
load are not as reliable as the incremental changes in strain for each incre- 
mental change in load. Normally, the tightening of the end clamps on the speci- 
mens resulted in an initial strain of some significance. 

The assistance of Mr. W. H. Taylor in carrying out the testing program is 
acknowledged here. 


Stresses 

Based upon the available mechanical properties as previously mentioned, 
the stresses were calculated from the stress-strain relation and the trans- 
formation relations. The tangential stress component ctq obtained from the re- 
corded strains are plotted in figure 4. 


FINITE ELEMENT ANALYSIS 


The finite- element method used here is a two-dimensional hybrid approach. 
The variational principle used is that of minimum complementary energy with 
the interelement stress continuity enforced by means of the Lagrange multipliers. 
The elements used are shown in figure 5. 

The formulation of the problem at this stage follows a rather standard 
fashion as this method is typically applied to many stress analysis problems. 

The stress function polynomial used in the computer program is 

i|f = ax^ + bx^y + cxy^ + dy^ (2) 


which results in the following stresses ’ 


<j = 2 cx + 6 dy 

XX ■' 


a = 6 ax + 2 by (3) 

yy 


T = -2 bx - 2 cy 
xy 

It can be easily verified that these stress components automatically satisfy 
the equilibrium equations in the absence of body forces. 
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Arranged in matrix form, equations (3) become 
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By Cauchy’s relation = a^.n^ , one has 
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Following a somewhat standard fashion, the circumferential stress around 
the hole is obtained and plotted also in figure 4 for the comparison purpose. 


DISCUSSION OF RESULTS 


The results presented in this paper represent an attempt to understand 
and predict the stress concentration around a circular hole in an angle-ply 
laminate. As shown in figure 4, the circumferential stresses^ based upon the 
finite- element method and the one computed from the recorded strain data, are 
plotted together for con^arison purpose. These two curves cross each other at 
a few places, but the discrepancy at some places is up to 357c This degree of 
deviation is not hoped for, but it is tolerable. Similar experience indicates 
that such a difference is by all means possible. 

The stress concentration factor at b/r = 1 and 0 = + 90° is about 5.8 
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which is considerably higher than the classical factor 3 for an infinite, 
isotropic plate. Therefore, special attention must be paid for the local 
stress concentration around such a circular hole. One possible reason for 
having such a high stress concentration factor is that a number of fibers were 
cut at the location of the hole. This weakens the ability of the fiber ele- 
ments for transmitting the stresses. From an intuitive point of view, if 
the location of the hole is known in advance, then rerouting the fibers 
around the hole may cut down the high stress concentration factor. 
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TABLE 1.- TEST SPECIMEN DATA 


Station 

Thickness, mm (in.) for specimen — 

1 1 

2 

3 

4 

0° 

■fiSSSni 

KE^H 


HPH 

45° 


wtg^gM 

■oES^fli 

■jHQQj^jH 


o 

o 


moQ^mi 

imogn 

HEESnE 

0.5842 

(0.0230) 

135° 

- 


■Q^IQIIIII 

■qKBH 

■||Q^|g|||[| 

0.5842 

(0.0230) 

180° 

HELfibH 

mmBM 

■oE^nB 

0.5969 

(0.0235) 

-135° 

WcEyiiM 

{■RjjQQIH 

HQgC^|H 

■SEgE» 

0.5969 

(0.0235) 

-90° 

■qEESH 

BuEl^n 

wg^gsgM 

McE^M 

0.5969 

(0.0235) 

-45° 

■EEBIH 

flCESSH 

0.6223 

(0.0245) 

0.6223 

(0.0245) 


TABLE 2.- TEST RESULTS 


Applied 

load, 

N (lb) 

Remote 
stress, 
kN/m^ (ksi) 

Strain recorded, mm/mm or in./in. , for specimen — 

HI 

2 

3 

4 

448.8 (100) 

6895 (1) 

1520 X 10 6 

-990 X 10 6 

390 X 10 6 

108 X 10^ 

889.6 (200) 

13790 (2) 

2610 X 10 6 


600 X 10 6 

250 X 10 6 

1334.4 (300) 

20685 (3) 

3730 X 10 6 

-2690 X 10 6 

830 X 10 6 

398 X 106 

1779.2 (400) 

27580 (4) 

4920 X 10 6 

-3540 X 106 

1020 X 10 6 

551 X 10 6 

2224 (500) 

34475 (4) 

6170 X 10 6 

-4600 X 106 

1230 X 10 6 

710 X 106 
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